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Abstract 

We study meson loop corrections to the baryon magnetic moments starting 
from a QCD-based quark model derived earlier in a quenched approximation 
to QCD. The model reproduces the standard quark model with extra cor- 
rections for the binding of the quarks. The loop corrections are necessary to 
remove the quenching. Our calculations use heavy baryon perturbation the- 
ory with chiral baryon-meson couplings and a form factor characterizing the 
structure of baryons as composite particles. The form factor reflects soft wave 
function effects with characteristic momenta ~ 400 MeV, well below the usual 
chiral cutoff of ~ 1 GeV. The resulting model involves only three parameters, 
the quark moments [i u and [i s and a parameter A that sets the momentum 
scale in the wavef unctions. We find that this approach substantially improves 
the agreement between the theoretical and experimental values of the octet 
baryon magnetic moments, with an average difference between the theoretical 
and experimental moments of 0.05/Ujv- An extension to the decuplet states 
using the same input predicts a moment of 1.97 fijy for the Or hyperon, in 
excellent agreement with the measured moment of 2.02 ± 0.05 
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I. INTRODUCTION 



The simple, nonrelativistic quark model (QM) gives a qualitatively good description of 
the baryon moments. Under the assumption that each baryon is composed of three valence 
or constituent quarks in a state with all internal orbital angular momenta equal to zero, the 
moments are given by expectation values of the spin moment operators 

M QA/ = £ w*> (i-i) 

q 

leading to the standard expressions 

^ M =\{^ u -^\..., ^=2^, (1-2) 

where the effective quark moments \i u = —2fid and fi s can be treated as free parameters in 
attempting to fit the data. The predicted pattern of the signs for the moments agrees with 
observation. A least-squares fit to seven measured octet moments (the transition moment 
/xeoa is left to be predicted) gives a root-mean-square deviation of theory from experiment of 
0.12 /xjv, about 10% of the average magnitudes of the moments. Agreement at this level can 
be regarded as an outstanding success of the quark model, but the deviations also give a very 
sensitive test of baryon structure. There is presently no completely successful first-principles 
theory of the moments. 

In this paper, we approach the moment problem dynamically using a QCD-based quark 
model |]] with meson loop corrections. The model has only three parameters, namely the 
quark moments \i u = —2 fid and fi s , and a parameter A that characterizes the momentum 
scale of meson-baryon interactions, with the particles regarded as composite. We find that 
this approach reduces the average deviation of the calculated moments from experiment to 
0.05/ijv, a substantial improvement, and that the loop corrections are small compared to the 
leading QM terms, suggesting reasonable convergence for the loop expansion. 

A different approach to the moment problem using chiral perturbation theory (ChPT) 
with loop corrections has been studied recently by a number of authors. Excellent fits to 
the data can be obtained for either conventional [T2],[3],[|,|5| or a modified U chiral counting. 
However, as pointed out elsewhere 0, the general fits are essentially independent of the dy- 
namical input, and are best thought of as giving parametrizations of the data. In particular, 
the introduction of the counterterms needed to eliminate divergences in the loop integrals 
results in the appearance of five new chiral couplings at the one-loop level in addition 

to the two tree-level || couplings. The seven well-measured octet moments can be fitted 
exactly using these seven parameters, whether or not the calculated loop corrections are in- 
cluded J7| . The theory is only weakly predictive with respect to the remaining quantity, the 
E°A transition moment, which is not known precisely. Furthermore, the loop corrections in 
ChPT are nearly as large as the tree-level terms, and convergence of the expansion is at best 
slow, when the divergences are regularized using dimensional regularization. This problem 
is not unexpected, since chiral symmetry is known to be badly broken for the baryons. The 
remaining predictions |2|,0,f§] are sum rules that can be motivated by large N c expansions 
or expansions in the symmetry breaking mass parameter m s in ChPT, but do not depend 
explicitly on the dynamics. 
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Our intent here is to present a dynamical approach to the baryon moments in which 
we emphasize the composite nature of the baryons and use a description based on the QM 
rather than the chiral picture. We begin in Sec. II by reviewing briefly the derivation of the 
QM, with corrections, from QCD using a Wilson-loop approach |l[|. The derivation involves 
the suppression of internal quark loops ("quenching"), an approximation that is likely to 
account for the deficiencies of the resulting model. The addition of quark loops necessary to 
remove the quenching introduces meson loops when the theory is viewed at the hadron level 
|J, and the effect of these loops on the baryon moments must be considered. We develop 
our approach to loop effects in Sec. Ill, where we define the chiral couplings we will use 
for low momentum meson-baryon and electromagnetic interactions in a calculation based 
on heavy-baryon perturbation theory PJlOH. We also introduce the form factor needed to 
account for the extended structure of the baryons and mesons viewed as composite particles. 
We present our final expression for the octet moments, the results of our numerical analysis, 
and our conclusions in Sec. IV. Finally, we give the detailed results needed to evaluate the 
formal expressions for the moments in two appendices. 



II. BARYON MOMENTS IN A QCD-BASED QUARK MODEL 

A. Background 

In previous work [|lj, we analyzed the QM for the baryon moments in the context of QCD. 
Our approach was based on the work of Brambilla et al. ||1 1|| , who derived the interaction 



potential and wave equation for the valence quarks in a baryon from QCD using a Wilson- 
line construction. Their basic idea was to construct a Green's function for the propagation 
of a gauge-invariant combination of quarks joined by path ordered Wilson-line factors 

U = Pexpigj A g ■ dx, (2.1) 

where A g is the color gauge field. The Wilson lines sweep out a three-sheeted world sheet of 
the form shown in Fig.l as the quarks move from their initial to their final configurations. 

By making an expansion in powers of l/m q using the Foldy-Wouthuysen approximation, 
and considering only forward propagation of the quarks in time, Brambilla et al. were able 
to derive a Hamiltonian and Schrodinger equation for the quarks, with an interaction which 
involves an average over the gauge field. That average was performed using the minimal 
surface approximation in which fluctuations in the world sheet are ignored, and the geometry 
is chosen to minimize the total area of the world sheet subject to the motion of the quarks. 
The short-distance QCD interactions were taken into account explicitly. Finally, the kinetic 



terms could be resummed. The result of this construction is an effective Hamiltonian fTI 
to be used in a semirelativistic Schrodinger equation Hty = Efy, 

2a s 



H = J2ypf + m i+ a ( r i + r 2 + r s) - o , . 

la, . 1 r , N «<= N a s 



-si • (ri x pi) + — -jsi • [(r i2 x Pi)— + (r i3 x Pl ) s 
2 1 a s 2 1 a s 

r s i ' r i2 x P2 - 77 3- s i' r i3 x P3H • (2-2) 
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Here r»j = x, — Xj is the separation of quarks i and j, is the distance of quark i from 
point at which the sum r% + r 2 + r 3 is minimized, and Pi and Sj are the momentum and spin 
operators for quark i. The parameter a is a "string tension" which specifies the strength of 
the long range confining interaction, and ct s is the strong coupling. The terms hidden in the 
ellipsis include tensor and spin-spin interactions which will not play a role in the analysis 
of corrections to the moment operator, and the terms that result from permutations of the 
particle labels. The full Hamiltonian is given in [[□]]. This Hamiltonian, including the terms 



omitted here, gives a good description of the baryon spectrum as shown by Carlson, Kogut, 
and Pandharipande |T2| and by Capstick and Isgur [T3|], who proposed it on the basis of 



reasonable physical arguments, but did not give formal derivations from QCD. 

The presence of the quark momenta in the Thomas-type spin-dependent interactions 
in Eq. ( |2.2| ) suggests that new contributions to the magnetic moment operator could arise 
in a complete theory through the minimal substitution 

Pi -> Pi - e i A em (a; i ), (2.3) 

with A em (xj) the electromagnetic vector potential associated with an external magnetic field. 
We have therefore redone the calculation of Brambilla et al. with the gauge interaction 
extended to include A em . After reorganizing the calculation to keep the presence of A em 
explicit throughout, and then expanding to first order in e q A em in the baryon rest frame 
with 

A em — — B x x 9 , B = constant, (2-4) 

we can identify the modified magnetic moment operator through the relation 

AH = -n-B. (2.5) 

The new moment operator, = ^t(Q M ) + A/x^ M , involves the leading corrections to the 
quark- model operator associated with the binding interactions. A/x can, in fact, be read off 
from the terms in Eq. (|2.2|) which depend on both the quark spins and momenta by making 
the minimal substitution in Eq. ( |2.3| ). For example, the term which depends on sj • r 12 x pi 
gives an extra contribution 



6^f Xl X (Sl X ri2) ^ ( 2 - 6 ) 



to [i v There are also possible orbital contributions to the moments because the Hamil- 
tonian mixes states with nonzero orbital angular momenta with the ground state. These 
contributions proved to be negligible 

The moment of a baryon b is now given by 

^ = £0 9 + A/i^K,,), = $ M + £ (2.7) 
q i 

where the sum is over the quarks in the baryon and we have quantized along B, taken along 
the z axis. The spin expectation values are to be calculated in the baryon ground state. 

Note that the correction A// to the moment of quark q depends on the baryon b in which 
it appears. The final baryon moments depart from the quark model pattern only when the 
ratios Afi b /fi q differ in different baryons. 
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B. Results for the octet baryons 



The general result for the moment operator given above can be simplified considerably 
for the L = ground state baryons. The absence of any orbital angular momentum allows us 
to integrate immediately over angles. Furthermore, two quarks always have the same mass 
in each octet or decuplet baryon, so appear symmetrically in the spatial part of a flavor- 
independent wave function. We will label these quarks 1 and 2, let mi = m 2 = m, and take 
3 as the odd-mass quark if there is one. With this labeling, the spatial matrix elements of 
the new operators can be reduced to a small set, and the correction terms become 



QM,b 



Mi 



/i 2 



/'3 



e + e' + 1_ + es_\ _ E + £' 
2m ei Vm m 3 / 2m 



2m e 2 Vm m 3 



E' 



2m 



e 

— + 

m 3 e 3 



1 ei + e 2 ( j 



m 



E 

m 3 



(2.8) 



where the e's and E's are ground state radial matrix elements, 

2a s r i2 ■ xi 

\ n „3 h 



( 



and 



2a s 


ri2 • xi 


9 


^12 


2a s 


r 23 ' x 2 


9 


/yi3 

' 23 


2a s 


r 31 ' x 3 


9 


7"31 


= - 


r 12 ' x i 


6 


T\1 


a 


T23 ■ X 2 


- - 
6 


^23 


a 


r 31 ' x 3 


= 6 


T31 



), = ( 



9 r 12 
2a s r 13 • xi 
9 



^31 



2a s r 32 ■ x 3 v 

\ n _3 /ft ' 



9 



r3 
23 



(2.9) 



)ft 
)ft 
)ft 



o- , r 31 • xi 
6 r 3i ! 



^" / r 32 • X 3 

6 r 23 



(2.10) 



In writing these results, we have made the approximation r\ + r 2 + r 3 « |( r i2 + r 23 + 
r 31 ), known to be reasonably accurate for the ground state baryons fll2H , and used the 
corresponding Thomas spin interaction. 

We have evaluated the the radial matrix elements above using Gaussian wave functions 
obtained in a variational calculation of the ground state energies for the Hamiltonian in Eq. 
( p.2|) . The results are given in Table | for a s = 0.39 and a = 0.18 GeV , values taken from 
fits to the baryon spectrum [[j3,|13[| using the same Hamiltonian. A refitting of the octet 
moments using fi u and fi s as adjustable parameters in the QM contribution to the complete 
expression in Eq. fl2.7| ) gives only a slight improvement in the results, with an incorrect 
pattern in the corrections relative to those needed. 
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III. LOOP CORRECTIONS TO THE MOMENTS 



Given the failure of the binding corrections to eliminate the deficiencies in the QCD- 
based quark model, we have reexamined the approximations used in its derivation A 
key element in the Wilson-loop construction was the use of the quenched approximation 
in which all internal quark loops are omitted. As discussed elsewhere JIJ, it appears that 
this approximation is the one most likely to account for the difficulty in reproducing the 
measured moments. In particular, the introduction of quark loops allows meson loops to 
appear, and these are known to affect magnetic moments. The first step in the improvement 
of the model is therefore the introduction of meson loop corrections to the QM moments. 
The relevant one-loop Feynman diagrams at the hadronic level are shown in Figs. |2| and || 
Since the ground state L = octet and decuplet baryons only differ in their internal spin 
configurations in the simple QM and the octet-decuplet mass differences are small, we must 
include both sets of baryon states in the calculation to get a consistent theory. However, we 
include only the pseudoscalar mesons. 

The diagrams in Fig. || are independent of the input magnetic moments of the baryons, 
but modify the final moments. In contrast, the diagrams in Fig. [3] involve the octet and 
decuplet moments and the octet-decuplet transition moment directly. We will specify these 
in terms of the QM. 

We will calculate the contributions of the diagrams in Figs. |2] and |^ using heavy baryon 
perturbation theory (HBPT) in which the baryon masses are assumed to be large compared 
to the typical scale set by the internal momenta. It will be essential in this respect to take 
the extended, composite structure of the baryons into account, since this extended structure 
naturally limits the the momenta that can be absorbed by the baryon as a whole. The 
resulting momenta are small on the average, well below usual chiral cutoff at ~ 1 GeV, and 
it is therefore reasonable to use the lowest order chiral couplings to describe the resulting 
the low momentum, long distance interactions of the mesons and baryons.S We will discuss 
these ingredients separately in the following subsections. 



A. Heavy baryon chiral couplings 

Heavy baryon perturbation theory was developed in Ref. || and extended to the chiral 
context in Ref. [|10| . It has been used to study a number of hadronic processes at momentum 



transfers much less than 1 GeV. The key ideas in HBPT involve the replacement of the 
momentum of a nearly on-shell baryon by its on-shell momentum itibV^ plus a small 



x The difficulties with ChPT in the present context result from the treatment of the baryons 
and mesons as point particles in HBChPT. This leads to divergences in loop integrals. Different 
regularization schemes lead to different results for the integrals, for example, in dimensional [14] 
and momentum regularization, with the ambiguities being lumped into the new couplings that 
must be introduced along with the loop diagrams. The present theory has a natural cutoff at 
fairly low momentum imposed by the extended structure of the hadrons in the QM, and no new 
couplings appear. 
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additional momentum k^, p = ttibv + k, and the replacement of the baryon field operator 
B(x) by an operator B v (x) constructed to remove the free momentum dependence in the 
Dirac equation, 



BJx) = e imB ^ x »B( 



x 



(3.1) 



In these expressions v M is the on-shell four velocity of the baryon, and it is assumed that 
k ■ v <C w,b- The perturbation expansion then involves modified Feynman rules and an 
expansion in powers of k/rriB- We will work to leading order in this expansion. 

The chiral Lagrangian for the modified baryon fields depends on the usual matrix of 
baryon fields, with B replaced by B v , and on the pseudoscalar pion octet normalized as 



1 

V2 



/ Zl + JL 
7T~ 



7T + 




K+ 


K° 


u 


K° 

2 V 
V6 



\ 



(3.2) 



This couples to the baryon fields at low momenta through the vector and axial vector currents 
defined by 



v» = r 2 ( w - w) + • • ■ , a, = r l 9^ + ■ ■ 



(3.3) 



where / ~ 93 MeV is the meson decay constant. We will retain, as indicated above, only 
leading term in the derivative expansion. The lowest order chiral Lagrangian for octet and 
decuplet baryons is then 

C v =iTrB v (v V) B V + 2D Tr B v S% {A, t , B v } + 2 F Tr B v S% [A^ B v ] 

- i T£ (v ■ V) T v , + 5 T* T VIX + C (T» A, B v + B v A^ Tf) 



+ 2Hf» S vu A" T V n + Tr d^cj) + 



(3.4) 



where 5 is the decuplet-octet mass difference, and = <9 M + [V^, ] is the covariant chiral 
derivative. B v is now the matrix of octet baryon fields, and the Rarita-Schwinger fields T£ 
pQj] represent the decuplet baryons. D, F, C, and 7i are the strong interaction coupling 



constants. The spin operator S% is defined in Ref. ||10 



The electromagnetic interactions of the mesons, and the convection current interactions 
of the baryons, are introduced into the Lagrangian by making the substitutions 



(3.5) 



where A^ is the photon field. The baryon moment couplings needed in the graphs in Fig. 
do not appear above, and require a separate treatment. 



B. Quark model moment couplings 

As described above, our starting point for calculating the baryon moments is the QCD- 
based quark model. These moments play the role of tree-level interactions in the loop 



7 



expansion, and also appear in the electromagnetic interactions of the internal baryons in 
Fig. ^|. As we noted in 0, the QM baryon moments can be written in terms of the general 
SU(3) symmetry breaking operator in HBChPT given by 

C SB = C + Ci (3.6) 

where 

C = -^L D TiB v F^{Q,B v } +fi F TrB v F tlv a^[Q,B v }) , (3.7) 

4777, n v ' 

is the leading order moment operator in chiral perturbation theory parametrized by /i£> and 
Hp ||, and 

Ci = -r—Fuu f ciTr B v MQa fJ,v B v + c 2 Tr B v Qa^ u B v M + c 3 Tr B v a^B v MQ 

+ c 4 Tr B v Ma^B v Q + c 5 Tr B v a" u B v Tr MQ ) (3.8) 

contains the new couplings C\, . . . , C5 introduced along with the counterterms that are nec- 
essary at one loop [FJ,[3],[|,|5]|| • Here Q = diag(2/3, —1/3, —1/3) is the quark charge matrix, 
and M. = diag(0, 0, 1) is proportional to the mass matrix used to introduce SU(3) breaking 
through the strange-quark mass m s . The addition to £33 of one further coupling which is 
second order in m s ,i 

ed 

£2 = 1 TrB v MQa tlu F^ByM (3.9) 

gives a complete basis for the description of the octet moments 0. 

The particular choice //p = 2//d/3 = fi u and c\ = — 5A, C2 = 0, C3 = —A, C4 = 0, and 
C5 = A, d = for the parameters above, with A = (2/x s + /i«)/2, gives us the QM moments 

4 1 

3^ s 3^ u ' 

14 . 

+ > ( 3 - 10 ) 

3 

Since the one loop corrections involve both intermediate octet and decuplet baryon states, 
we also need the Lagrangians for the decuplet magnetic moment couplings and the octet- 
decuplet magnetic moment transitions for the QM. These are given respectively by 

£(3/2) = -i^-T^ m Q)T v J M F^ , (3.11) 



..(QM) _ 
H'p 


3 

2^ ' 


..(QM) _ 
PS+ ~~ 


4 


1 


„(QM) 


..(QM) _ 
H'n 


H"U 1 


(QM) 
/^ E o — 






Ms- 


..(QM) _ 

AA - 




,,(QM) _ 


2 

-3^- 


1 

"3"" 


„(QM) 



and 



£(«0 = -i^ F ^^ Q\Bi m S%T» klm + fc.c) , (3.12) 



where i, j, fc, I, and m are SU(3) flavor indices and Q = diag(/x u , — /i u /2,// s ). 



2 Since AiQ = QA1 = — gA^f, the rearrangements of the factors AfQ and Af in Eq. ( |3.9| ) give no 
new contributions to £2, and the form given is the only new invariant. 
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C. Meson wave function effects: The form factor 



The baryons in the QM are composite states, and can absorb only limited recoil momen- 
tum while remaining in their ground states. This must be taken into account in a dynamical 
model. In the absence of a detailed theory of the interactions of composite mesons and 
baryons, we will model the wave function effects using a form factor at each meson-baryon 
vertex. In keeping with the heavy baryon picture, we will define the form factor in the rest 
frame of the baryon where it can depend only on k 2 , the square of the three momentum of 
the meson. A form factor of this type automatically repects crossing since a change k — > — k 
corresponding to a shift of a meson between the initial and final states does not change the 
form factor. 

We have chosen for simplicity to use a one-parameter form factor 

m = < 313 > 

normalized at chiral limit, k = 0. The parameter A characterizes the natural momentum 
scale, expected to be much below 1 GeV. With the introduction of this form factor, all the 
diagrams in Figs. ^ and |3] are finite and no arbitrary subtractions need to be introduced into 
the theory. 

The evaluation of momentum integrals in the presence of the form factor is straightfor- 
ward, but involves some new elements. First, the form factor (|3.13| ) is rewritten in covariant 
form as 

F <M = g _ ( _ A , ■ < 3 - 14 > 

The factors in the denominator of a one-loop integrand with integration variable k can be 
combined using the Feynman parametrization formula to obtain an expression of the form 

k 2 + a(k-v) 2 + P(k-v)+-f . (3.15) 

Here a, ft, and 7 are quantities independent of the loop momentum k, and the photon 
momentum q has been set equal to zero in the denominators. A change of the variable of 
integration to 

k' = k + [±y/l + a - 1] v(k ■ v) (3.16) 
then brings Eq.( [3.15| ) to the standard form 



k' 2 ±^=(k'-v) +1 , (3.17) 
V 1 + « 

and the loop integral is easily evaluated. Note that the Jacobian of the transformation of 
variables in Eq. (|3.16|) is 1/ y/1 + a. 
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IV. BARYON MAGNETIC MOMENTS 



A. Expressions for the octet baryon moments 

We can now give our expressions for the baryon magnetic moments including the loop 
corrections from the diagrams shown in Figs. |2| and [| In units of nuclear magnetons, the 
moment of baryon i is 

A^^ +M f /2) +/xf /2) , (4.1) 
where the leading term /4 includes the QM moments plus the corrections A^f M obtained 

(1/2) (3/2) 

in the QCD-based QM, while /i- and /i- are the contributions from the one-loop graphs 
with intermediate octet and decuplet baryon states, respectively. We find that 

/xf =«< + , (4.2) 

(i/2) = m N A 4 ( x) 

* xkic 24vr/MA + m x )3 ft 

+ E ^ Y ^(ll iX) -2X ( i X) a l )L (m x ,X) , (4.3) 



and 



(3/2) _ V- m N ~ ~ {X ) 



+ E ^Ar2 tf?* ) -2>l* ) a i )L 1 (m x ,8,\) + j^ ) L 2 (m x ,8,X) , (4.4) 



X=n,K 

^9^2 n L 

where = ■ The coupling coefficients , A- X ^ + X$ are identical to those 

in Ref. ||], and will not be given here.! The remaining coefficients r )i' X \ ^l^-, and j^ X ^ 
depend explicitly on \i u and /i s , and have not appeared elsewhere. We list those coefficients 
in Appendix A. 

To connect the various terms to the loop graphs in Figs. ^| and || we note that (3\ X \ 
fti i ll i ll i ll i ^\ X \ an d are the coefficients of the graphs 2a, 2b, 3a, 3b, 3c 
(or 3d), 3e and 3f, respectively. 

The expressions for the functions L (mx, A), F(mx, <5, A), Li(mx, 8, A), and L 2 (mx, <5, A) 
are given in the Appendix B. These functions result from the loop integrations, and depend 
on the meson masses, the decuplet-octet mass difference 5, and the natural wave function 
cutoff A. We emphasize that the loop integrations are all finite, and that the wave function 
parameter A sets the natural momentum scale in graphs that are divergent for point particles. 



3 The tadpole graph 2b in Ref. [Q] is absent in our model. Since the chiral symmetry is broken in 
the QM, we do not include chiral corrections to the baryon moment operators as in that reference. 
The elementary BB(p(pA^ vertex connected with spin operator therefore does not exist. 
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B. Fits to the data 



We have used the foregoing expressions to fit the experimental baryon magnetic moments. 
It is worth emphasizing that all the parameters that appear except for fi u , fi s , and A are 
known independently. We evaluated the corrections A^f M given in Eqs. ( p. 71) and Q2.8|) 
from the QCD-based QM using the matrix elements given in Table H and and quark masses 
m u = m d = 0.343 GeV and m s = 0.539 GeV chosen to give the best fit to the octet magnetic 
moments in the naive QM. The strong interaction couplings F, D, and C were chosen to 
satisfy the SU(6) relations F = 2D/3, C = —2D expected in the QM, with the values 
F = 0.5, D = 0.75, and C = —1.5 chosen so that F + D « \qa/9v\ — 1-26. The decuplet 
coupling TC does not appear. The decuplet-octet mass difference was taken as 5 = 250 MeV; 
the results are fairly insensitive to this choice. Finally, we used the values f n = 93 MeV and 
/* = /„ = 1.2/,, @. 

The results of an equal weight least square fit to the seven well-measured octet moments 
using the three free parameters fi u , [i s , and A are summarized infj We find root-mean-square 
deviation of the predicted values from the observation of 0.055 fijy. This is a substantial 
improvement on the naive quark model which gives an average deviation of 0.12 fi^. The 
transition moments /is^A was not included in the fit, but was left as a prediction. The result 
obtained, /i£°A = 1.559 /x^, differs from the experimental value /iso A = 1.610 ± 0.08 /ijv, 
by 0.05 /xjv; a value within the experimental uncertainty and one that does not affect the 
overall fit. 

A detailed breakdown of the contributions of the loop integrals to the fitted magnetic 



moments of the octet baryons is given in the [TTll The results in this table shows that the 



loop contributions are small in comparison to the leading QM contributions, suggesting 
reasonably rapid convergence of the loop expansion. This is in marked contrast to the 
results obtained in HBChPT, where the loop contributions calculated using dimensional 
regularization are comparable in size to the leading terms in the chiral expansion J7|. It 
is also clear from the table that the binding corrections Afif M found in the QCD-based 
derivation of the QM |] are important. It is not possible to obtain as good a fit to the data 
as given in the Table [TI] when these are omitted. 

A further interesting point involves the importance of the decuplet intermediate states 
for the octet moments. It is easy to check that the contributions from the graphs involving 
decuplet states, that is, the sum of the contributions from the graphs 2b, 3b, 3c, 3d, and 
3f, are substantial. For most of the baryons, those contributions are larger than those from 
the graphs which involve only the intermediate octet states. This result is insensitive to 
the value used for the octet-decuplet mass difference. We conclude that the decuplet must 
be regarded as a set of light baryon states as in the QM, and not, as in as a set of 
heavy states. The present method for obtaining the loop corrections to the QM moments 
should therefore apply to the decuplet moments as well. That calculation has been carried 
out elsewhere jllj using the same values of fi u , and A as obtained here, and gives a 



prediction = —1.97 /xjv in striking agreement with the measured value —2.02 ±0.05 /xjv, 
and much better than the prediction of the naive QM, /j,q- = —1.74 fi^. 
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C. Conclusions and prospects 



In this paper, we have considered the one-loop corrections to the octet baryon magnetic 
moments in the QCD-based QM. It is necessary to include these as a first step in getting 
away from the quenching of internal quark loops used in the derivation of the QM. 

Our approach to the calculation is based on HBPT. We use the derivative couplings 
for the low-momentum interactions of the pseudoscalar mesons with the baryons favored 
by ChPT, and a form factor to characterize the composite structure of baryons. The loop 
diagrams are all finite, and no counterterms need to be introduced to absorb divergences. 
We are, in fact, making a dynamical calculation of the extra couplings or counterterms 
encountered in HBChPT in the sense that our expressions for the eight octet moments can 
be parametrized exactly |7|] in terms of the eight chiral couplings defined in Eqs. ( pj.7|) , ( p\8| ) 
and Q. 

The results from our fits to the octet baryon moments using the QCD QM with loop 
corrections are very good, with an average deviation of the theoretical moments from theory 
of 0.05 fix, significantly better than the QM results. The contribution from each loop 
graph is small compared to the leading terms, suggesting convergence of the loop expansion. 
The parameter A which sets the momentum scale in the meson-baryon interaction or wave 
function is about 400 MeV, a value consistent with the expectations deduced from the 
observed transverse momentum distributions in pion production. This value is closer to the 
kaon mass than to the pion mass. As a result, the wave function effects suppress the short- 
distance contributions from kaon and ^-loops, but affect the more reliable long-distance part 
of the pion loop contributions relatively little. We conclude that it is crucial to take the 
effects of compositeness into account if one is to have a controllable perturbation theory for 
hadrons at low momentum. 

A question which arises at this point concerns the extent to which the theory can be 
further improved by the inclusion of higher-mass mesons and baryons. The contributions of 
the ground state vector mesons would be expected to be as important as the contributions 
from the pseudoscalar mesons except for the suppression of high mass intermediate states 
by the momentum cutoff imposed by the wave functions. The contributions of higher mass 
baryons in the intermediate states are also suppressed. However, there are low-mass mul- 
tiparticle intermediate states such as those with one baryon and two pions that could be 
important. A possible approach to the estimation of their effects is through the use of the 



sideways or mass dispersion relations proved by Bincer [17|. 

We remark, finally, that we believe that this work demonstrates the importance of getting 
beyond the quenched approximation in lattice QCD if one is to understand the finer details 
of hadron structure from first principles. 
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APPENDIX A: THE COUPLING COEFFICIENTS 



In this appendix, the coupling coefficients are explained and presented explicitly. For 
simplicity, the superscript (X) is suppressed. Our fa, and fa are identical, respectively, to 
the coefficients fa and in Ref. Q. The sum of our coefficients Aj and Aj is equal to the 
coefficient A,- defined in 0, 



Aj + A» — Aj. 



(Al) 



Since the couplings associated with graphs which involve only octet intermediate states are 
independent of those for graphs which involve decuplet intermediate states, it is easy to 
separate A, and Aj from the combined coefficient given in [[J. 

The values of the coupling factors jj evaluated from the graphs in Fig. 3a are 



ip 

In 
1(tt 

7a 

1(tt 

7s+ 

1(tt 

7 E o 

1(tt 

7e- 

1(tt 

7 H o 

1(tt 

7h- 

1(tt 

7aeo 



-(D + F) 2 ^ u , 
2 

^[(3DF-5F 2 )fi u + (2F 2 



D 2 )fi s 



g[-2F 2 /in 



(2F 2 - D 2 )/i s 



-[(-3DF + F 2 )fi u + {2F 2 - D 2 )fi s 

-2(F-D) 2 fz s , 

^F-D) 2 ^ u -8fi s ), 

1=D(4F-D)fi u , 



(A2) 



for the pion loops, 



7 



1(K) 
V 

l(K) 



7 

l(K) 

7a 

1(K) 

7s+ 

1(K) 

7 s o 

1(K) 

7 E - 

1(K) 
Too 



(i? -D)(D- 3F)/2 U + (^D 2 - 2DF - F 2 )fi s 



2F(F - D)pi u 



-D 2 - 2DF - F 2 )fi s 



-^(D 2 + \2DF + 9F 2 K - ~(D - 3F) 2 fi s 
to y 

- 7 -(D 2 - 2 lDF + F 2 K - i(D + F) V, 

_I (jD 2 _ 4jDF + f2) _!(£> + F) 2 /i s , 
o 3 

_ —DF + F 2 )^ - -(D + F) 2 /i„ 
o5 3 

-2D(D + FK + i(7J 2 + 6-DF - 3F 2 )/i s , 



(A3) 
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7 ig ) = ^ 7 =(3D 2 + 4DF-9F 2 K 



7^ (??) = -t(3f-l>)V 



for the kaon loops, and 

1 

4 

7 ; ( "> = ^(3F - 2?) v, 

7i (,?) = -\d^ S) 

l 1 ^ = - 2 -D\^ u -^ s ), (A4) 

7^ = ^(2^ + ^), 

7 i ( o ??) = ^P + 3F) 2 (^-4 / , s ), 

7H W = -^P + 3F) 2 (/i« + 8/i s ), 

7ae° = 7y|^ A*« 

for the rj loops. 

The coefficients 7/ evaluated from the graphs 3b are given, up to a factor of — 5C 2 /2, by 

~1M 16 ~1W 1 /r , /I \ 4 

7/ ; = - -qI*u, 7s+ = -^( 5 ^« + 4 ^)> 7 a o = -gA*«, 

7* W = gV«, 7$° = -^(^ + 2^), 7= ( - ) = + 4/i s ), (A5) 

7a W = -\(^u + 2fi 8 ), 7s- } = ^(A*u - 4^ s ), 7aeo = 
for the pion loops, by 

2uA ^l (K) = — (23u„. + 4//.V ^li K) = --(3/i u + 14/i s ), 

9 

1 



Ip 


-^(3/i„ + 2/i s ), 


7^ = 


- — (23/^ + 4/^), 


7 S <> = 


- 

In 


1 . 

-{Hu -2fi 8 ), 




-^(13/i u + 8/i s ), 


ll [K) = 


~1(K) 

7a = 


1 . 




4 . 

— {5fi u - 2/i s ), 


7ae<> — 



(A6) 



for the kaon loops, and by 



ip 


= 0, 


7#> = 


In 


= 0, 


7#> = 


-1(7,) 

7a 


= 0, 





3^ 



2 2 

■-(2/i„ + /i s ), ^ = -g(Vu + 2/i S ), 



gvr-- r~»n , a - g 

g(/i«-/is), 7aeo = 0> 
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1 (/i u -4/x s ), (A7) 



for the r] loops. 

The coefficients 7 2 evaluated from the graphs 3c (or 3d) are given, up to a factor of 2C, 

by 

7 P 2W = ~(D + F)n u , fff = ~l((3D + 5FK - 8F/i s ), 7^ = \{F - D)^ 

7* W = l(D + FK, 7^ = ~F{Hu - 4/i s ), lP = \{D - F)(5pt u - 2^), (A8) 

7i W = - 4/x.), 7s W = ^((3I> + F)fi u + SF^s), 7 S = + 6^K 



for the pion loops, by 



f P {K) = -\[D^ u + {F-D)^ 
fW = l[(D + F) P L u + 2(D-F) P L s ], 
^ = ^(3F-D)(ji u -A f ji a ), 
7$° = ^[(F- 5D)pt u + 2(D + 

nS C) = ~(D + F){^-A^) 1 (A9) 
^ = ^[(2D-F)^ + (D + F)^ 3 ], 
7$° = ^[(£> + 3F)/i M - 2(£> + F)fi s ], 
^1 {K) = -^[F f i u +(D + F) f i s ], 



7^ = -^(2£> + 3FK 



for the kaon loops, and by 



1? ] = 0, lP = ~D{^ - /i s ), 7^ = 9 (D + 3F)(/i tt - /, s ), 

2 „ , . o ( 2 



7^=0, f^ = --D(v u -4» s ), 7? ) = -^P + 3F)(/i u + 2/i s ), (A10) 
7^=0, 7^ = ^(^ + 2/i s ), f^ = -^=D^ 
for the 77 loops. 

APPENDIX B: THE EXPRESSIONS FOR L , F, L u AND L 2 

Let us introduce the following notation 

\/m 2 — a 2 [7r/2 — arctan (a/ V"^ 2 — a 2 )] m> a 

F Q (m,a) = I (Bl) 

— \J a 2 — m 2 In [(a + V^ 2 — m 2 )/m] m < a 
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where a is an arbitrary parameter. Hereafter, the subscript (X) is ignored. The function 
L (m, A) obtained from the Feynman integral for the graph 3a (or 3e) is then given by 



L (m, A) 



A 



^(A 2 + 2m 2 ) + Am 2 F (m, A) 
3 \ z — m z 



(B2) 



(A 2 — m 2 ) 2 

The function F(m, 5, A) obtained when calculating a Feynman integral for the graph lb is 



7rF(m, 5, A) 



X s 



+ 
+ 



3(A 2 -m 2 + 5 2 ) 2 
X5 



N(m,5,X) + 



5A 2 + 2m 2 „ A 2 + 2m 2 



-5 + 



(X 2 -m 2 ) 2 (X 2 -m 2 + 5 2 ) 
3m 2 



A 2 — m 2 (A 2 — m 2 ) 2 
3(2A 2 + 3m 2 )(A 2 - m 2 ) - 2(A 2 - Qm 2 )5 2 



A 2 — m 



:5> 



F (m,A) , 



(B3) 



where 

N(m,5,\) 



(A 2 - m 2 + 5 2 ) 



7rA(A 2 + 3m 2 -3S 2 ) -2(3A 2 + m 2 - <5 2 ) F (m, 5) . (B4) 



Similarly, the functions Li(m, 5, A), and L2(m, 5, A) arise from the Feynman integrals for the 
graphs 3b (or 3f) and 3c (or 3d), respectively. We have 



Li(m, 6, A) 



2\ 4 



+ 



3(A 2 -m 2 + 5 2 ) 2 l 

A 

(A 2 -m 2 )(A 2 -m 2 + <5 2 ) 



(A 2 + 2m 2 - 25 2 ) - 5N(m, 5, A) - ^ 9 A - < 

\ z — m z 



2 X 2 



3m 2 (A 2 - m 2 ) - 6X 2 S 



2A 2 - 3m 2 
A 2 — m 2 



5 4 



F {m, A) 



(B5) 



and 



L 2 (m, 5, A) 



2\ 4 



+ 



3(A 2 -m 2 + 5 2 ) 2 

7T<5 



2(m 2 - S 2 ) F (m, S) - 7rm E 



2(A + m) 
A 



/a 2 a a 2n A + 2m r2 
(A + Am + 4m 2 ) - - <r 



A + m 



A 2(A 2 - m 2 + £ 2 
A 2 — m 2 



(A 2 — m 2 ) 2 
When 5 = 0, Eq. (IBJ) gives 



3m 2 (A 2 - m 2 ) - (2A 2 - 3m 2 )5 2 F (m, A) 



F(m,0, A) 



A 1 



3(A + m) 3 

and it follows from Eqs. (|B5|) and (|B"6|) that 

Li(m, 0, A) = L 2 {m, 0, A) = 2L (m,A) 



(B6) 



(B7) 



2A 4 



3(A 2 -m 2 ) 2 



[(A 2 + 2m 2 ) + -|^F (m,A) 
\ z — m z 



(B8) 
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FIGURES 



FIG. 1. World sheet picture for the structure of a baryon in the Wilson-loop approach. 
FIG. 2. Diagrams with couplings independent of the baryon moments. These diagrams lead 

1/2 

to the non-analytic m s corrections to the baryon magnetic moments in the conventional ChPT. 
The dashed lines denote the mesons, the single and double solid lines denote octet and decuplet 
baryons, respectively. A heavy dot with a meson line represents a form factor F(k,v) (Eq. (|3.14| )), 
where k is the meson momentum. 

FIG. 3. Diagrams with couplings that that depend on the tree-level baryon moments. These 
diagrams lead to non-analytic m s In m s corrections to the baryon magnetic moments in the con- 
ventional ChPT. 
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TABLES 



TABLE I. 


The values in GeV of the matrix elements 


e and E defined 


in Eqs. (Ej9 


) and ((2.10) 


The matrix elements were 


evaluated for a s 


= 0.39, a = 


0.18GeV 2 , m u = 


= m d = 0.343 GeV, and 


m s = 0.539 GeV. 












Baryon 


e 


e' 


e 


£ 


£' 


£ 


N 


0.014 


0.014 


0.014 


0.062 


0.062 


0.062 


E 


0.014 


0.018 


0.012 


0.061 


0.069 


0.046 




0.016 


0.012 


0.018 


0.053 


0.046 


0.068 


n 


0.016 


0.016 


0.016 


0.052 


0.052 


0.052 



TABLE II. The magnetic moments from the naive quark model (QM) and the QCD-based QM 
with loop corrections. The average deviations from experimental values are about 0.10 and 0.05 for 
the QM and the QCD-based QM with loop corrections, respectively. All moments and deviations 
are given in units of hn- 



M-B 


QM 


QCD QM + loops 


Experiment 


P 


2.728 


2.720 


2.793 ± 0.000 


n 


-1.818 


-1.946 


-1.913 ±0.000 


£+ 


2.618 


2.519 


2.458 ± 0.010 


E~ 


-1.019 


-1.110 


-1.160 ±0.025 


E° 


0.800 


0.705 




A 


-0.580 


-0.608 


-0.613 ± 0.004 


H° 


-1.380 


-1.316 


-1.250 ±0.014 




-0.470 


-0.582 


-0.651 ±0.003 


S°A 


1.575 


1.559 


±1.610 ±0.08 




1.818 


2.083 




Ms 


-0.580 


-0.656 




A (in MeV) 




407 
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TABLE III. Detailed breakdown of the contributions of the loop integrals to the fitted magnetic 
moments of the octet baryons (in jUjv)- Those contributions are evaluated at F = 0.5, D = 0.75, 
C = —1.5. As shown in Table I, a best fit is obtained at [x u = 2.083, fx s = —0.656 and the natural 
cutoff jj, = 407 MeV. stands for a deviation from the experimently measured value. 





Mu> Ms 


Afj,i 


ml /2iN) 


In mi N ^ 


ml /2iA) 


In 


Loops 


HB 




p 


3.124 


-0.500 


0.416 


-0.729 


0.057 


0.353 


0.096 


2.720 


-0.072 


n 


-2.083 


0.427 


-0.381 


0.409 


-0.067 


-0.253 


-0.291 


-1.946 


-0.033 


£+ 


2.995 


-0.428 


0.273 


-0.536 


-0.002 


0.217 


-0.048 


2.519 


0.061 


E" 


-1.170 


0.139 


-0.217 


0.188 


0.021 


-0.072 


-0.080 


-1.110 


0.050 


S° 


0.913 


-0.145 


0.028 


-0.174 


0.010 


0.073 


-0.064 


0.705 




A 


-0.656 


0.085 


-0.028 


0.073 


-0.010 


-0.073 


-0.037 


-0.608 


0.005 


E° 


-1.569 


0.271 


-0.043 


0.128 


-0.043 


-0.059 


-0.017 


-1.316 


-0.066 




-0.528 


-0.056 


-0.048 


0.061 


0.033 


-0.044 


0.002 


-0.582 


0.069 


E°A 


1.804 


-0.378 


0.229 


-0.268 


0.058 


0.114 


0.134 


1.559 


-0.051 
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